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SUMMARY 


Chaplygin introduced the hodograph method in the theory 
of compressible fluid flows and developed a method for oxr ^ 
constructing, stream' functions of such flows.., T-hi s method - ,^ 
which has teen -extensively used in investigation of compress- 
ible fluid flows, is limited in certain respects-:-- The ex- 
pression for the etreajo function obtained in this manner, can 
represent only- certain typep of flow patterns, : Ip, general, 
flow patterns obtained in this way cannot represent the whole 
flow around an obstacle,' but only a par-t- of such a -flow, and 
therefore several expressions are. needed in order- to obtain 
the whole flow.- On the other hand, in many -instances, it is 
important to have a ¥ .single expression r epr c sept.ing the whole 
flow, . • ■ • • - 

Recently Von Karman and Tsien constructed more general 
tyoes of stream functions, but only by replacing the true 
pressure density relation. ■■ - - ... 

p = ap k 

by the linear pr e s sur e- s-p'e cif i c volume relation 

p = A + — 

.-. ■ ■ P 

(A, k, and U are constants), so that their method is essen- 
tially limited to flows the maximum Mach number of which is 
not too large. ' 

In a companion report the author' derived a new 'formula 
for stream functions based on the true pressure density'r e-" 
lat i on 
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p * ap k 


It is not subject to the limitations of the Chaplygin method. 

In the present report this formula is employed to con- 
struct two-dimensional subsonic c ompr e s 8 i bl-e fluid flowB 
around a body similar in shape to a given symmetric obstacle. 

The methods described in the report are illustrated by 
numerical examples. J 


INTRODUCTION 


In'a companion report (reference l) the author derived a 
f ormula 1 transforming an arbitrary analytic function of’ a com- 
plex variable, g(z), into a function \|/(x-,y) which satis- 
fies the differential equation for the stream function' of a: 
two-dimensional potential gas flow. (It - is assumed that- the 
pressure of the fluid Is a function of its density.) * 

Thi s report describes the application of this formula 'to 
the determination of a flow around an obstacle. First-, the 
actual’ c omoutat i on of the stream function generated by -an 
analytic function is. described in detail. Then various meth- 
ods for choosing the function which yields a flow around a 
given obstacle (or at. least around a shape similar to the 
given' one) are discussed. <" •. 

The procedure is illustrated by completely- carrying o.ut 
the construction of the compressible flow generated by the 
function 


g(Z) = 


= 1 ( 3 

2Vyr~; 




2e 


2e ‘ 


r ) 


In the case of an incompressible fluid this function 
wou-ld lead to the' circulation-free flow aro.und a circular 
cylinder. In the case of a .compressible -fluid a flow around 
a cylinder of somewhat distorted cross section .is obtained. 

For the mathematical- justification of the formulas used, 
the reader is referred to reference 1. 
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This investigation, conducted.. at -the Brown University, - 
was sponsored "by and conducted with the financial assistance 
of the National Advisory Committee for Aeronautics. • 

I take the opportunity to express my gratitude to Mr. 
Leonard Greenstone for his assistance in the preparation of 
the present paper. ' 


NOTATION 


a speed of sound 

a 0 speed of sound at a stagnation point 

f(z) analytic function of the complex variable z 

g(Z) analytic function of the complex variable Z 

g ( n+i ) ( 2 ) nth iterated integral of g(z) (See (16).) 

H a function of a real variable occurring In the formula 
for the stream function of a compressible flow (See 
sec. 2.) 

k exponent in the pressure-density relation 

M Mach number (M = v/a) 

( n ) 

Q; function of a real variable occurring in the formula 
for the stream function^ of a c ompr e s sdble flow 

( n ) * 

T n imaginary part of g considered as a function of M 

and 6 

v speed (magnitude of the velocity vector) 
v 1 ,v 3 components of the velocity vector 

w = cp + i\|/ complex potential 

z = x + iy complex variable in the physical plane 

Z = A - i8 complex variable in the logarithmic plane 

angle between the velocity vector and the x-raxi s 


6 
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\ for an i nc ompr e s si tie flow: logarithm of the speed.; f-or 

a compressible flow: a given function of threr local 

Mach number (defined in sec'. 2) (The X of this re- 
port ia not to be confused with A of reference 1, 
to which it is similar, but differs from it by a con- 
stant.) ... 

cp velocity potential 

cp n real part of (See (20).) 

\j/ stream function 

\|/ imaginary part of (See (20).) 


1. Construction of a Flow by Means of— Analytic Functions of a 
Complex Variable (Incompressible Fluid) 

The method used in this paper consists of generating 
two-dimensional compressible flaws by means of analytic func- 
tions of a complex variable. The method is best explained by 
first considering the case of an i nc orapr e se i bl e fluid. 

There exists a very simple method of constructing pat- 
terns of a two-dimensional, 'steady irr otrat i onal flow of- an in- 
compressible perfect fluid "by means of analytic functions of 
a complox variable. Let 

W = f(z) = cp ( x , y ) + i \j/( x , y ) , z = x + iy (l) 

\|/ satisfy the Laplace 

sfvk + , o 

dx 3 By 3 

where cp can be int erpret ed _ as' 'the velocity potential of an 
incompressible flow; \j/ is the stream function, and 

\|/(x,y) = constant 
is the eauation of a streamline. 


be such a function, 
equat i on , 


Then cp and 


2 . ,3 . 

a co + o cp _ 

8x 3 By 3 


and 


For instance, 
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x(x s * y 2 + l) + ^ (x 3 + y 3 ‘ - i) 
2(x s + y 2 ) 2 ( x 2 + y 3 ) 

* .» • i 


( 2 ) 


represents a circulat i on-free flow around a circular cylinder, 
for Im(w) is constant for . jz| = 1. (See fig. 1.) 

Unfortunately, thi s, ,si.jnple method cannot be extended to 
the case of a compressible fluid. . . 


Now, let and- v 2 he. .the velocity components of the 

incompr e s sihle flow . 


' ; ■ - v x = dcp/Sx, v a = dcp/dy , 

) 

Let v "be the speed and 8 the angle between the velocity 
vector and the x-axis. Then 

; " ' : - iVa = ve“ i& = (3) 

Thus v x - iv 2 is also an analytic function of z. By vir- 
tue- of (3) to each point (x,y) of the physical plane, there 
corresponds a point (v lf ~v 2 ) of the plane the Cartesian 
coordinates of which are v x , -v 2 (hodograph plane). In 
this way ah. image of- the flow in the (v lf -v 2 )-plane is ob- 
tained. This Ima^e is called the hodograph of the flow. It 
may happen that at two different points of the (x,y)-plane 
(.physical plane) v x and y 2 have the same value. Then, 
to these two points there corresponds the same point in the 
hodograph plane. If It is' desirable to have a one-to-one 
correspondence between the 'flovand it g -hodograph ,- the latter 
must be interpreted aB a "Itiemann surface." (See reference 
2, pp . 27', "59~$0, etc.) ' 

1 ' * 

J 

Tor the flow ( 2 ) 



The .hodograph of the flow consists 
circle ( v - I J + v 2 3 ' = i. ■ Each 

v 2 r ' 4 

corresponds to two 'points in the 


of, all points within the 
point of the hodograph 
(See fig. 2.5 


■z-plane . 
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In the hodograph plane v and -0 -are polar coordi- 
nate Sv ; 

The "logarithmic" plane is defined as the plane in 
which A , 

A * log v 

and -6 are the Car t e s i an , c o or d i nat e s . The complex variable 
in the logarithmic plane is 



The transformation (4) maps the hodograph of the flow into 
its image in the logarithmic plane. Note that the mapping of 
the hodograph plane into the logarithmic plane is independent 
of the f 1 ow . 

For instance, the image of the flow ( 2 ) in the logarith- 
mic plane is the (doubly covered) domain bounded by the curve 

A = log 2 + log cos 0 , -^ 8^7 (•&) 

(See fi.gs. 3 and 5.) - ’■ ' 

From the .basic properties of analytic functions it fol- 
lows that the complex potential w. ' also, can he consicLere.d as 
an analytic function of Z, 

’ w = f ( 2 ) = g ( Z ) ; . ' t 6 ) 

Fop instance-, fcr "the flow ( 2 ) ........ 

Z = iog (l ~ \y, that is, z = (l'~ 2e*)~^ ( 7 ) 

and therefore 

g(Z) = [(1 - 2e Z j^ + (1 - 2e Z )"^ (8) 

The streamlines in the physical plane (,-z-plane) corre- 
spond to the lines ’in the logarithmic plane along which 

C ■ ’’ v T 0 (A ,0 ) = 1m g(z ) = \j/ ( 9 ) 


i s constant . 
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F-or' tthe flow (2) 


T 0 (X-,6) -= lm 

_ 1 

2j2 


g(Z) 


- (1 - 2e X cos G) + .(1 - 4e X cos 0'- + 4e 2X )^ 
- X [l - 1 1 


(l - 4e^ cos 0 + 4e 2 ^)^-* 


( 10 ) 


(The values of T 0 are given in table III and plotted, in 
fig. 5.) In this .example the streamlines are the same in 
the' two sheets of -the Eiepann surface. 


Now, let an ar bitr ary ..apaly 1 1 c function g(z) he given. 
It always is possible to interpret this function as the com- 
plex potential (in the logarithmic plane) of an incompressible 
flow in the physical plane (z-plane). The flow in the ? r plane 
can be determined as follows. Since w(z) = g(z) t and . 


then 


dz dg -Z . 
dZ ~ dZ 6 


or 



Integration along a streamline (in the Z-rplane) yields a 
streamline in'the z-plane. In this case the foregoing inte- 
gral can be written in the form „ " 



In order that the flow be physically possible, it is necessary 
that the streamlines do not intersect. 
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Application of this procedure to function (4) leads back 
to the flow ( 2) . 

By means of the foregoing method of associating with 
every analytic function the Z-plane (logarithmic plane), a 
flow in the z-plane (physical plane) can be extended to the 
case of a compressible flow. 


2. Construction of Flows by Means of Analytic Functions of 

a Complex Variable (Compressible Fluid) 

Consider 'a' two-diraensi onal subsonic potential flow of a 
compressible fluid (in the S-plane), It is assumed that the 
density p and the pressure p are connected by the relation 

• p = A + arp^ ( 11 ) 


where A, o', and k r are constants. There exists a stream 
function \| / (x,y) such that 


ip,. 

P dy 


. 1 9i|/ 

p Bx 


( 12 ) 


are the components of the velocity vector in the x and y 
directions, respectively-.”, 

Let M be the local Mach number 


M = M ( v ) 



(13) 


where a * 0 

the flow, 
positive 


is the speed of sound at rest, v the speed of 

0 the angle between the velocity vector and the 
x-axis. Set 



^ 1 - . ( 1 - M s ) ^ 

1 1 + h(l - M 3 )^ 

1 ' 
h 

j 

^1 + (.1 - M 3 )% j 

|l - h(l - M s )^j 

_ 


2A = log 


(14) 
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where ' 



k - 


k + 


_ 1 * 
. 1 _ 


k > 1 


(If k < 1 the definition, of X must he changed. (See ref 
er ehce 1 . ) . 


The logarithmic plane is defined as the plane with 
Cartesian coordinates X, 0. The complex variable in this 
plane is 

Z s* X - i 0 


In the Z~plane \}/ , which may be considered as a function 
of X and 0, satisfies the equation 

lb ( \J/ ) = I + N M 5 §!i + I M B 0 (15) 

0 4 Vaxa aea.y ax Bza"z Vsz az/ 

where — 

N = -Qc + 1)M 4 

8(1 - M 2 ) 3/s 


This is a linear eauation and its treatment can be facili- 
tated by using results from the highly developed theory of 
linear partial .differential 'equations. -Let g(z) be an ar- 
bitrary analytic function of Z. Using this function, _a 
solution of (15) pan be obtained as follows: 


Set 


g (o) \z) = g(Z) 


g (n+l) (z) = r g (n) (.Z)dZ, n = 0, 1, 


( 16 ) 


There exist fixed fun'ct ions -of a real variable | 


H(£). Q Cn> (l ) . n « 1. 3. 


such that the funotion 
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H(2X - 2a) 


Im g(z) + \ - 2a) I m g (l> (z) 

(p 

.(3) 


+ 1 C^ 8 \2A - 2a)lm g v * ; (z) + . . . 


+ llsil Q, C n ) ( _ 2 a)lm g (n)( Z ) + 


2 Sn ni 


'(17) 


(a a real constant) satisfies (15). (The representation (17) 
is an approximation to the solution. The exact formula 
proved in reference 1 is 


\l/ (^“.9^ = H( 2 X - 2 a) 1 i 9 
\- a O / m 


m — >a> L 


Im g( z ) 


Y <^J n ' > (2\ - 2 a)lm g (n) (z) 


n=>i 


3n f “ 
2 n ! 


' (18) 


where a is an arbitrary non-negative constant 


lim Q, 
m — >■ od 


(n) rt (n) 


m 


In the. example under consideration i,t is replaced by (17), 
since the Mach number is comparatively small and the 

do not differ consid.er-ably from o/ n ^ . If, however, the 
local Mach number approaches the value 1, it is necessary to 
use the exact formula. (Another form of formula (17) is 
given in appendix I.)) The functions H and Q,( n ) depend 
only upon k and v. They have been defined by a recurrence 
formula and computed for n = 1, 2, 3. (See reference 1, 
sec. 12.) Some of them are tabulated ‘in tables la and lb and 
plotted in diagrams la and Ih . . 

• f 

Once \j/ has been determined as a function of v/a Q and 

8, the t r an.s i t i on- to the physical plane (z-plane) is given 
by the formulas . 
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x 


y 




c. o s 6 M 3 \\iq . 


pv 

sin 8 M 2 \J i 


dv 


> (19) 


fl d 


dv 


pv 


where the integration is performed along a streamline (i.e., 
along a line \|/ = constant) and the subscripts 8, v denote 
partial differentiation. ~ " 


In this way an arbitrary number of streamlines in the 
physical plane (z-plane) can be drawn. It should be noticed 
that along each streamline the value of v/a D i b known and 

therefore also the values of the local Mach number, ‘the pres- 
sure, and the density. 


The next section contains a discussion of the actual 
computation of the function \J/ and of the streamlines in the 
physical plane. " ■ . 


Two additional remarks may be made; 

1. The series (17) converges rather slowly for values of 
v in the neighborhood of v = a (i.e., M = I).. - This is 
not due to a singularity of ty, since it is known that \|/ 

is regular in the domain concerned. Therefore it is possible 
to replace (l?) by other expansions which converge more rap- 
idly. (See reference 1, sec. I?-) The author also. will treat 
'this question in a future publi cat i oiijr ... . .... 

2. It Is known that the theory becomes especially simple 
if k = -1. (The pr e s sur e- spe c i f i c volume relation then be- 
comes linear.) Therefore it is interesting to -know how 
(generated by a given analytia function g(z)) changes if k 
varies in the interval (-1,2). For this reason the computa- 
tions are carried out for two values of k, k= -0.5, and 

k = 1,4 i'h'the example given in the following section. 

3. Actual Computation of the Stream Function, Using the 
Generating Qp.et.atOr in the Form (17) 

The actual computation of a compressible flow given in 
the logarithmic plane by an analytic function g(z) consists 
of the following steps: , 
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I. Computation of the function g( n )(z) 

II. Construction of the flow in the logarithmic plane 

III. Transition from the logarithmic plane to the physi- 
cal plane 

These steps are discussed in detail in the present sec- 
tion. The procedure is illustrated by constructing the com- 
pressible flow obtained from the function 


g(Z) 



(1 - 2e Z )“^ 




(1 - 2e Z )^ 1 


( 8 ) 


(Por the case of an incompressible flow this function yields 
the flow around a circular cylinder; see the first section.) 

St en I The functions g^ n ^(Z) have been defined in 
the preceding section. In general, it will be impossible to 
evaluate the integrals in closed form 


g ( n+1 ) ( Z ) = j g (n) (Z)dZ 

X+i6 „° 

cp (n) (X,6) + i* (n) (X,0) 


(dX - id 6) 


X , 6 “ ‘ • * • X , & 

(cp (n) dX '+ \|/ (n) d0) t.i ■ 


(-jqp^ n ^d0 + \|/ n ^dX) 


( 20 ) 


where 


9n + M*i = e 


(n) 


however, it "is nossible to use numerical integration. Since 
g (n) (z) is an analytic function of a complex, var iable , the 


value of the integral 


g( n ) (z)dZ 


is independent ofhthe 


% 

path of integration. .A very convenient 'path consists of the 
interval [(0,6), (X,~6)]' of the negative real axis, followed 

by the Interval [(X,o), (X ,9 )] along a line parallel to the 
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0 axis (provided this path lies within the image of the flow 
in the 1 ogar ithm.ic /plane . " If this, dobs not take plape, L an_ 
obvious adjustment, must be made.) "Dividing each of these in- 
tervals into s subintervals AX k , A6 k gives approximately 



s • , 

- --'i -V'- cp k < n >A 0 k . (21) 

k^i 


I'n many instances it is more convenient to e xpe^n.d g ( Z ) 
into an infinite series and then integrate term by term. 
However, as a rule, function g(Z) has singularities, and 
several., series developments are heeded in order to cover the 
domain in which the function is to be considered, For 'the 
case .of function (8) integral / gdZ can be comuuted in 
closed form. Then 



Z ' 

g( Z)dZ 


(1 



log 



1 + (l 



0 , *■ — -! . 

' - Z - i (^2_|_rr^ ( 22 ) 


and 

T 1 (X,0) = Im 


1_ 

72 


M 

J s(Z)AZ 




-(l - 2e^'cos 0 ) + (l - 4e^ cos 8 + 4e^) s ~ \ ~ ( — 2~^) + 2 


. ... »>. i. . ... 

■ -iv i 3 


f 2^ e^ Bin 8 + [-(l - 2e^ cos 8) (l - 4e^ cos 8 + 4e^)^] 

^ 2^('e^ cos 0 - l) — . £(l - 2e^ cos 9) + (l - 4e^ cos 6 + 4e^)^3 


( 23 ) 


The values of T 0 (X,6), (10), and of T 1 (X, 8) for various 
values of' (X;,8), are given* in t.ahl.e II*.. ... - r. 

■ • “ - . 1 - ■*- * *• t 

The functions g, T 0 (X,8)> and T x (X,8) can be repre- 
sented. by. the; following series-. Her r 0.691 < X < 0: 
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i/2 


' Z /2 3 -Z/2 -3Z/2 

e - e - — — e — — ■ — e 

4 32 128 


5Z/2 




. 9 - 7Z /3 _ 


2048 


T 


o " *S 2 


* /Z 008 (I)-- I e ’ X/S cos (I) - Js e ’ 3> ' /2 e0B 

- t. r 5?>/2 cob (a) _ _m_ e - ,x / 2 cos 

128 ■ \ 2 / 2048 


(¥) 


° S (jf) + • ’ *_ 


>(24) 


=a /2 


32 


' ^ / 2 " , 


3 2 / -A/2 

e cos 

1) - \) 

f + — X— ( e cos 

-3A/2 

cos 

<£)- 

i> + _Z_ x i(',- 5 V 2 
/ 128 5 V 


+ -IS- x 2 (" e - ?X /3 

2048 7 \ 


cos 


jr A < .-0.691s 
= - 2 


(i 2 ) ■ 0 + • • ■. 


t , 1 2Z . 3Z , 15 4 Z 

1 + — e + e + r 


+ e 4Z + I e 5Z + 
8 2 


S« = -2 


2A 


3A 


2 • 


e sin 28 + e sin 36 + iS e" ' sin 46 


4A 


8 


.7 5A 

+ — e " sin 50 + . 
2 


H25) 


T 1 * *- 2 


& + ^ e 2 ^ s i n ( 26 ) + i e 3A sin(36) + || e 4 * sln(46) 


+ — e 5X sin (50) '. . 


+ 1.14 


J 
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St et> II . - The coftstr’aot ion of the. flow in the logarith- 
mic plane can he oonv.en.ien.tly performed on. specially scaled 
graph paper. One "(Car t"e sian-)- coordinated axi s is taken as 
the 8-axis, . the other as'the M-axis. In addition to the 
values of M s'cales showing, the .corresponding values of 
A(M) and v/a 0 should he indicated. Such scales may he 
preoared once and for all, for each given value of k. (See 
(ll) and (13).) In diagrams' la and Ih. these scales are 
drawn for k = -0.5 and k = 1.4. 


The values found for : -Im 
and the lines 



are entered on this paper 


Im g^ (A 


i8) 



onst ant 


‘ (26) 


f - ’/ l 

are drawn. In this way there is obtained by graphical inter- 
polation the values of the functions 


, Un(M l 0) = 



X(M) 



(27) 


The function ty(M,0) is determined by formula. "(17)^ .. 
In applying this formula- if is necessary to choose a. t de’f’4. nit e 
value for the arbitrary 1 constant a. In general-, it is., pos- 
sible to choose a so' that the terms in the series (17) be- 
come .small for large, value, .of M. . In the example under con- 
sideration a is set equal to 0.1.' 


In order. to obtain the function \|/(M,0 : ) it is necessary 
to evaluate the products * ’ . 


| q (l) (M)T 1 (M,0 : ), ■ | q (2) '(M)T 2 (M,8), ... (28) 

- - r • * 

The values of the functions Cj/^(M) can be tabulated once 
and for all (for a given ■tfalue of ■ ‘k)- r f \The above products 
can be evaluated graphically,' by means "of '.a simple . aomogr.am, 
shown in figure 4. 


Adding a finite number of the terms in the series (16) 
and multiplying by H(M) (a function tabulated once and for 
all) (see tables. Ia, lb) gives an approximate expr e sM on for 
'|/(M,0). Finally, the lines \j/(M,B) = constant (stream 
lines) are drawn. 


For function (8) \|/(M,0) has ;beah. .,c omput ed. for k = -0.5 

and k = 1.4. For k - -0.5 the series ha'B been replaced by 



16 


NACA TN No.- '973 


\MM,0) = \|/( M , 6 ; - 0.5) 

= Im js( 2 X - 0 . 2 ) T 0 (M,e) 

+ 1 O' 1 ^ ( 2X - 0.2)Ti (M, e ! )l) , X = X(M) (29) 

2 ■' ‘ ' j J 

For k = 1.4 the series has "been replaced by 

\|/ Q - ( M , 0 ) = ' \|/( M , 0 ; 1.4) • ;■ 

. = Im (h( 2X - 0.2) T 0 (M,a) ,+ . Q, ^ ^ ( 2 X - 0 . 2 )T 1 (M, 6) 

+ £ <^ 2) (2X - 0.2)T a (M,e)l ] , X = X(M) (30) 

" ' 4 - • 11 J 

The streamlines (in the logarithmic plane) are shown in fig- 
ure 6 (k = -0.5) .and figure 7 ,(k = 1.4). 

i • , • 

Step III .- The flow is first transferred from the loga- 

rithmic plane to the hodograph plane.' This is done by using 
(for the hodograph plane) polar coordinate graph paper. The 
polar, coordinates in -the hodograph plane are v/a 0 and 0.’ 
Since for each plane the values of v/a Q and 9 are known, 
the transition presents no difficulties. ... 

The hodograph • streamlines of the flow generated by the 
function (8) are shown in figure 8 (for k = -0.5) and In 
figure 9 (for k = 1,4). 

In order to obtain' the streamlines .in' the physical plane 
it is necessary to compute the integrals (19) along the 
streamlines in the hodograph plane. 

tfhe partial derivatives of entering ‘In the inte- 
grands cafi -be 'expresses ag follows ■ 

\|/ v = R* 0 ^' Im g + g + £ R^ : Im g (l ^ 

Z • ■ . . * . ' 2 

- (s) (s) 

■',= +£ R ■ I m g +:. . . (3l) 

4 ' . 

I 

. ___ * 

(Se.e reference 1, sec.. 12.) 
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where 

R (o) _ H dX_ ^(i) = -h_ + -A. . . . 

dv v 2 dv 


(n) 


- (• 2n 2 )1 


2n-2 


(n - 1) ! 




2n(2n - l) g-^Cn) 
4n 


■^>( 32 ) 
dv 1 ■ ■ 


n = 2,3, . . . 


..J 


and 

M/ e / 


The. values of bA° \ . . . , for- k »= 1 . 4 are given in 

t able. I b . :• 


The formulas for the- derivatives of T 0 , T x for the 
case of the function (8) are given’ in ; ap£endlx ’I I . The Val- 
ues of the derivatives of - \f/ are gi'v’en r: in table III. 

The function p ■ which enters in the integrand also can 
be tabulated once and for all. 


H- 


Be g r 


+ 15. q 
8 


<*> He .*<=> 


+ •. 


(33 


The integration is to be performed graphically. When a 
number of streamlines in the physical plane are drawn, a 
rather complete picture of .the 1 flow is obtained’. In fact ,, 
for each point in the physical plane for which. the correspond- 
ing point in the hodograph plane is known, the value of. the 
Breed, pressure, local Mach, number, and so forth, also is 
known . 

For the. function (8) it turns out that, a part of the 
streamlines \|/ = 0 forms a closed contour. The streamline 
starts at -as and 'divides, it self info two. branches at the 
first, stagnation; point . The. two .par t s come together at the 
second stagnation point . Thps a flow around an o.bstacle t has 
been obtained. 1 - The' boundary of the" obstacle (for k = -0.5) 
i s shown in f i gur e- 10 . ’ ’ 1 


x For k = 1.4, see p. 21. 
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4. Determination of the Function g(z) 

Which Leads to a Flow of Desired Type 

It has "been shown that every analytic function g(Z) 
leads to a compressible flow in the z-plane (physical plane). 
The -type of the flow obtained thus depends upon the choice of 

g(z) . • • 

If it is desired to obtain a flow around a given obsta- 
cle, or at least around an obstgclo similar to a given one, 
it is well to proceed as follows: 

First, determine the incompressible flow around the given 
obstacle and obtain the complex potent-ial as a function of 
log v - 16 . (By using Theodorsen's method (reference 3) this 
problem can always be solved.) This function, g(Z), is 
transformed into a stream function of a c ompr e e.sible flow. 
Assume that the compressible flow in the a-plane is a flow 
around a closed hody. 1 

Let B be the original profile, B x the boundary of 
the obstacle of the compressible flow. Clearly, B 1 will 

be different from 3. For instance, if B is a circle, there 
is obtained (for certain values of the arbitrary constant a 
entering in the computations) the profile shown in figure 10. 
The profile distortion can be represented by associating to 
each point a of B a point z of 3 x and writing 

z 1 = z + d 

(z, z x , d are complex numbers). Now, the profile B* 
the. point s of whi.ch are given by. 

. z * = . z - d 

can be described as being distorted "in the opposite direc- 
tion.” If' oho were' to construct the incompressible flow 

^•Necessary and- sufficient conditions for a function 
\|/(v/a 0 , 6 ) to lead to a' flow around a closed body are de- 
rived in reference 1, sec. 14 . If the original profile is 
symmetric with respect to both the real and the imaginary 
axis, the compressible flow pattern will always flow around 
a closed profile. 
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around B* , . oomnute its -complex potential in the logarithmic 
plane, g*(Z), and use 't'his funct i on "tf o .-o bt.ain a compressi- 
ble flow, it might be expected that a 'flow Around a profile 
B 3 would be obtained which better approximates ; B, This 
procedure may be repeated. ■ • ■ * ■ • 

In general, the function g(z) will be of a very compli- 
cated form, even if the complex potential in the physical' 
plane is a simple expression. Therefore, it is convenient to 
know functions which can be expressed in a closed form -by 
means of elementary functions, and- w’hich are potentials’ of-/', 
incompressible flows around closed profiles (in the physical 
plane). Bach such, function leads to compressible flows of 
similar type, and the computation of these- flows is relatively 
simple. , = 

A function of the desired type depending upon three arbi- 
trary parameters will presently be constructed. By the follow- 
ing. successive, transformations, g as a function of Z is 
det 0 .rmi.ned'!, . . *• ' " • 

g = iZ x (34) 

Z 3 •= i— . 0, < A < oo , 0 < C < oo C3 5 ) 

... C + Z ^ . 


z 3 = J z 3 2 + J z 3 4 + 3 Z S S (36) 

_ ■ m + 'iZ_- , 

Z 4 ■= -i log - £ -t — -. .m real. . (3?) 

....... m - iZ 3 

5 )‘- . («) 

Z = log Z_ . ; (39 ) 

* 'll " ’ . . '*■ ■ • 

Here, A ; 1 C , ni 'are 'arbitrary constants. The right half 
plane of the " Zj- -plane - is .taken -into the domains indicated 
in f i gur es llb td lid, and into 'the domain indicated in f i g— 
ure 12. Now, consider a Riemann surface consisting of two 
sheet s 'spread over -the dOSiain of the Z B -plane (fig. 12 ) and 
possessing a branch -point at the point -cc- (corresponding to 
g = cb) . Furict 1 6ri ■ -'gCz g ) ■ 'maps -thi s Riemann surface into the 
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exterior of a slit along the' real 'axis (in .the g-plane). 
Therefore the domain in Z 5 -plane‘ can he interpreted as a 
hodograph of a ciroulat i on-free flow around a slender Symmet- 
ric body with a' sharp' trailing edge, 'and g' is the complex 
potential of the flow (provided the. constants are chosen L.n 
an appropriate way). 

The . Z-plane i.s the logarithmic plane for the above- 
mentioned flow' and therefore g(Z) can be t-ransf ormed. into, 
a stream function of a compressible flow by the method given 
in the preceding sections. • 

The computation will be -facilitated by the relatively 
simple form of g. . - . : 

. CONCLUDING EBMAEK5 

1 *. , ■ . . . «■ 

This paper employs the hodograph method introduced- into 
the theory of compressible fluids by Chaplygin (reference 4) 
in 1902. 

An essential feature of the hodograph method is that the 
principle of superposition ofL solutions (each of which repre- 
sents a stream function in the hodograph plane) holds. If 
ij/x(v,6) f \|/ 2 (v,0 ) , . . . are solutions of the differential 

equation (17) and if A 1 ,’ A 2 , . are arbitrary constants 

'then \|/(v,6) = E A n \i/ n ("v , 0 ) is also a solution of (17), One 

of the main problems of the theory consists in determining 
functions \|/(v,9) which yield in physical plane flows around 
closed bodies all streamlines s t ar ti ng * and ending at infinity 

Until the present time, no solutions of this kind have 
been given except that of Ton 'Karman (reference 5) and Tsien 
(reference 6), who replaced the' correct pressure density re- 
lation p = Cp 1 * 4 by p = A + 0" / p . This simplification 
yields a rough approximation and' ca'n be employed only when 
the maximum speed is sufficiently small. 

In the present caper a' solution which fulfills the* con- 
ditions described above and which is based ’on the true pres- 
-sur a density relation is derived' f o:r the first : t“ime 

’ * t * 1 

Two ' exampl e s ha-v-e teen 'considered; one corresponding to 
a : nearly' circular domain', the other; corresponding to a slen- 
deh obstacle.' In the 'first case the computation of the 
stream pattern has been carried out completely. 
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• The foregoing method may he applied to other profiles hy 
either of two alternative schemes. Another analytic function 
g(Z) may he chosen and the procedure employed in the present 
paper repeated. A combination of simple solutions, such as 
those-’ ''given hy Chaplygin (reference' 4}-, Bar's and Gelbart (ref- 
erence 7), or the author (reference 8, p. 23 and reference 9, 

p. 277) may he added to_either of the stream functions jfound 

in the’ 'afore-mentioned examples. In each instance stream 
functions may he found that fulfill the previously described 
c ondi t ions. 

, j • r 

The method can he -extended to partly supersonic flows. 

One pf the advantages of the foregoing method is that a large 
part of the computatio'n is independent of the function g. 
These - computations can he carried out once and for all and 
tabulated, which will greatly facilitate the application of 
the method. 

As has been seen i‘n the example Under consideration, the 
actual calculation of the (subsonic) flow corresponding to a 
given function g(z) does not entail any thePretical diffi- 
culties. • On -the o.ther hand, it does involve- a very consider- 
able amount of' numerical computation. It Is, consequently, 
expedient to use special computational devices, such as punch 
card 'machines, the differential analyzer, and- so forth, in 
order -to overcome these difficulties. The method described 
in the present section has, however, been- d eve 1 op'ed ‘ wi t h the 
assumption that no such special devices (except perhaps an 
ordinary computing ( mul t ipli cat i on’ and addition) m'acbin^ are 
available . .r . 

In the sequel to this report the author will 'discuss the 
use to which these devices may be put in computing compressi- 
ble fluid flows. 1 : • 


Brown University, 

Providence, R. I., May 15, 1944. 


- 1 T-h e flow in the physical plane. -for k 3-1.4" will be 
described on occasion as well as several aerodynamic conse- 
quences which follow from the present results." 
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APPENDIX I 

Other form of Formula (l?) 


It is often more convenient to use an operator which 
transforms analytic functions into stream functions and 
which differs slightly from (17). This operator has the 
form 



This operator can he obtained from (17) by setting 

g(z) = r f Z(1 - t 2 )') d»t / ( 1 - t 2 )^ 

■-L , 3 * 


or 


f(Z) = — 

TT 


r 1 TT /a 

/ Z sin * i£lgg. ^ H dd + g (0) 

J d(Zsin 2 d) 77 
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APPENDIX II 

Formulas for the Derivatives of T for the Function (8) 


3T 


39 


o _ 1 

2 


(j^ - cos 6) + - e^ cos 0 + e 2 ^^ 


-iC 1 - 


e ^ cos 


+ fi cos 6 + e 

\4 


8 + e 2 ^ 

4 


-i -4 


e^ sin 9 


2X\t 


— - cos 0 + e 


[-8 - •’ 

2X 7“ 3 /a 




3 T o _ 
3X ' = 


a T o e 2A 

■cot 8 -57— + — — 
o Q p 


— - cos 8 +e 


2X1-4 


.1(1- 

2 \4 


e cos 8 + e 


2AJ 


-4 



3T_i 
38 " 


-(I - eA 008 e ) + (4 - 


^ cos 6 + ,S XNili 


1 + 


1 ( 7 k r- e^ cos 8 + e 2 ^ ^ 

2 V4 


)-4- 


3Ti 

~3X 


To 


COB 6 
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Table la 

(Jc = -0.5j *.= 0.1) 

Functions used in computing the stream functions 
of a compressible flow. 


2A-2»4 

1 

M 

y/a 

o 

H 


-2.82 

0.265 

0.272 

1.000 

-0.0299 

-1.96 

0.403 

0.430 

1.000 

-0.0602 

-1.612 

0.473 

0.518 

1.002 

-0.0909 

-1.51 

0.497 

0.551 

1.003 

-0.1115 

-1.388 

0.528 

0.593 

1.004 

-0.1267 

-0.79 

0.690 

0.861 

1.015 

-0.3009 

-0.394 

0.820 

1.165 

1.046 

-0.7033 

-0.188 

0.901 

1.440 

1.111 

-1.5871 

■ -0.132 

0.921 

1.527 

1.145 

-2.2682 


Table lb. - The values of F, H, Q^, R^ ( for k = 1.4 » °< = 0.1): Functions used in computing 


the stream- function of a compressible floe. 


21-2® 

U 

▼/a 0 

F 

H 

qW 

Q U> 

a <3) 

q (4) 

rW 

rW 

r< 2 > 

h (3) 

— eo 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

.0000 


.0000 

.0000 

.0000 

-3.8772 

.1000 

.0999 

.0001 

1.0000 

.0000 

.0001 

-.0001 

.0002 

9.9600 

.0010 

-.0020 

.0000 

-2.5096 

.2000 

.1992 

.0011 

1.0002 

-.0009 

.0015 

-.0013 

.0008 

4.9197 

.0030 

-.0103 

.0148 

-1.7327 

.3000 

.2972 

.0064 

1.0014 

-.0056 

.0085 

-.0083 

.0064 

3.2141 

.0090 

-.0411 

• 

.0670 

-1.2071 

.4000 

.3938 

.0256 

1.0042 

,.0199 

•0342 

-.0395 

.0374 

2.3371 

.0235 

-.1203 

.2307 

-.8238 

.5000 

.4 879 

.0866 

l.ono 

' -.0574 

.1160 

-.1680 

.2116 

1.7947 

.0520 

-.3152 

.7600 

-.6706 

.5500 

.5341 

.1565 

1.0167 

-.0935 

.2101 

-.3496 

.5211 

1.5901 

.0755 

-.5082 

1.4060 

-.5364 

.6000 

.5795 

.2839 

1.0247 

-.1501 

.3823 

-.7494 

1.3527 

1.4150 

.1088 

-.8277 

2.6949 

-.4204 

.6500 

.6242 

.5245 

1.0359 

-.2401. 

.7089 

-1.6847 

3.7802 

1.2620 

.1566 

-1.3796 

5.4353 

—3203 

.7000 

.6680 

1.0060 

1.0515 

-.3870 

1.3663 

-4.0770 

11.7558 

1.1254 

.2275 

-2.3946 

11.8166 

-.2207 

.7500 

.7110 

2.0623 

1.0811 

-.6959 

2.8304 

-12.7121 

51.4676 





-.1615 

.8000 

.7532 

4.6583 

1.1049 

-1.0896 

6.4088 

-35-1357 

194-4562 

.8847 

.5248 

-9.0769 

81.6373 

-.1015 

.8500 

.7945 

12.5662 

1.1517 

-2.0186 

17.4338 

-I49.6938 

954.643 

.7717 

.8743 

-21.9464 

306.6745 

-.0535 

.9000 

.8349 

46.6378 

1.2275 

-4-3787 

65-3783 

-997.3960 

16432.652 

.6555 

L.6SQ4 

-71.6427 

1772.3621 

.0000 

1.0000 

.9329 

OO 

OO 

— OO 


- CO 


-.0000 

©O 

— ©O 

©O 
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Computation of the stream-function (in the logarithmic plane) of a 
compressible flow generated by the analytic function (8). 


-A 

0 

To ' 

T 1 

, , T 2 

t 3 

KK =- 0 . 5 ) 

vKS = 1 . 4 ) 

0.02 

0.05 

- 0.0376 

- 0.0020 

0.0000 

+ 0.0000 

- 0.0365 

- 0.0400 

0.02 

0.20 

+ 0.0002 

- 0.0402 

- 0.0008 

+ 0 . Q001 

+ 0.0250 

+ 0.0105 

0.02 

0.40 

+ 0.0968 

- 0.1550 

- 0.0049 

+ 0.0021 

+ 0.1950 

+ 0.0994 

0.02 

0.60 

+ 0.2051 

- 0.3318 

- 0.0160 

+ 0.0103 

+ 0.4163 

+ 0.0702 

0.02 

0.80 

+ 0.2975 

- 0.5588 

- 0.0413 

+ 0.0326 

+ 0.6518 

- 0.3030 

0.06 

0.08 

- 0.1159 

- 0.0044 

- 0.0003 

+ 0.0000 

- 0.1162 

- 0.1212 

' 0.06 

0.30 

- 0.0251 

- 0.0892 

- 0.0058 

+ 0.0006 

+ 0.0144 

- 0.0187 

0.06 

0.34 

- 0.0019 

- 0.1148 

- 0.0074 

+ 0.0008 

+ 0.0496 

+ 0.0084 

0.06 

0.40 

+ 0.0345 

- 0.1582 

- 0.0105 

+ 0.0016 

+ 0.1062 

+ 0.0453 

0.06 

0.50 

+ 0.0964 

- 0.2418 

- 0.0174 

+ 0.0041 

+ 0.2069 

+ 0.1008 

0.06 

0.70 

+ 0.2096 

- 0.4494 

- 0.0390 

+ 0.0167 

+ 0.4155 

+ 0.1387 

0.06 

0.90 

+ 0.2973 

.- 0.7016 

- 0.0755 

+ 0.0466 

+ 0.6181 

+ 0.0038 

0.06 

1.10 

+ 0.3440 

- 1.0154 

- 0.1340 

+ 0.1042 

+ 0.8066 

- 0.4074 

0.10 

0.15 

- 0.1791 

- 0.0154 

- 0.0019 

+ 0.0000 

- 0.1870 

- 0.1852 

0.10 

0.35 

- 0.0625 

- 0.1200 

- 0.0127 

+ 0.0005 

+ 0.0219 

- 0.0574 

0.10 

0.50 

+ 0.0414 

- 0.2446 

- 0.0246 

+ 0.0029 

+ 0.1258 

+ 0.0536 

0.10 

0.60 

+ 0.1073 

- 0.3442 

- 0.0363 

+ 0.0078 

+ 0.2266 

+ 0.1130 

0.10 

0.70 

+ 0.1664 

- 0.4564 

- 0.0505 

+ 0.0135 

+ 0.3248 

+ 0.1652 

0.10 

1.00 

+ 0.2975 

- 0.8540 

- 0.1387 

+ 0.0683 

+ 0.5936 

+ 0.1361 

0.20 

0.22 

- 0.3504 

- 0.0184 

- 0.0071 

- 0.0003 

- 0.3497 

- 0.3505 

0.20 

0.40 

- 0.1802 

- 0.1468 

- 0.0255 

- 0 . 0017 . 

- 0.1534 

- 0.1853 

0.20 

0.70 

- 0.0013 

- 0.5354 

- 0.0768 

+ 0.0044 

+ 0.1177 

+ 0.0098 

0.20 

0 . 80 . 

+ 0.1259 

- 0.5972 

- 0.0931 

+ 0.0105 

+ 0.2599 

+ 0.1310 

0.20 

1.00 

+ 0.2204 

- 0.8790 

- 0.1323 

+ 0.0313 

+ 0.4179 

+ 0.2138 

0.20 

1.10 

+ 0.2544 

- 1.0328 

- 0.1481 

+ 0.0536 

+ 0.4864 

+ 0.2322 

0.30 

0.30 

- 0.4688 

+ 0.0306 

+ 0.0162 

- 0.0050 

- 0.4689 

- 0.4692 

0.30 

0.75 

+ 0.0099 

- 0 . 5364 

- 0.1034 

- 0.0008 

- rO.0886 

+ 0.0164 

0.30 

0.85 

+ 0.0753 

- 0.6758 

- 0.2058 

+ 0.0060 

+ 0.1750 

+ 0.0743 

0.30 

0.95 

+ 0.1287 

- 0.8222 

- 0.2398 

+ 0.0250 

+ 0.2504 

+ 0.1235 

0.30 

1.10 

+ 0.1904 

- 1.0546 

- 0.2491 

+ 0.0325 

+ 0.3467 

+ 0.1888 

0.30 

1.20 

+ 0.2209 

- 1.2164 

- 0.3078 

+ 0.0428 

+ 0.4013 

+ 0.2166 

0.40 

0.35 

- 0.5617 

- 0.0258 



- 0.5616 

- 0.5577 

0.40 

0.60 

- 0.2034 

- 0.3274 



- 0.1682 

- 0.1543 

0.40 

0.85 

+ 0.0091 

- 0.6802 



+ 0.0843 

+ 0.1111 

0.40 

0.93 

+ 0.0568 

- 0.8016 



+ 0.1457 

+ 0.1770 

0.40 

1.05 

+ 0.1147 

- 0.9900 



+ 0.2247 

+ 0.2632 

0.40 

1.30 

+ 0.1941 

- 1.4062 



+ 0.3504 


0.60 

0.40 

- 0.6396 

+ 0.0278 



- 0.6360 

- 0.6376 

0.60 

0.45 

- 0.5273 

- 0.0422 



- 0.5245 

- 0.5241 

0.60 

0.60 

- 0.3070 

- 0 . 274E 



- 0.2892 

- 0.2864 

0.60 

0.75 

- 0.1572 

- 0.5112 



- 0.1239 

- 0.1189 

0.60 

0.95 

- 0.0208 

- 0.8364 



+ 0.0336 

+ 0,0419 

0.60 

1.00 

+ 0.0051 

- 0 . 9256 



+ 0.0664 

+ 0.0745 

0.60 

1.15 

+ 0.0683 

- 1.1744 



+ 0.1446 


0.80 

0.50 

- 0.3970 

- 0.0318 



- 0.3959 

- 0.3962 

0.80 

0.70 

- 0.2158 

- 0.3878 



- 0.2022 

- 0.2148 

0.80 

0.90 

- 0.0898 

- 0.7378 



- 0.0640 

- 0.0880 

0.80 

1.05 

- 0 . 0209 , 

- 1.0040 



+ 0.0143 

- 0.0169 

0.80 

1.10 

- 0.0017 

- 1.0934 



+ 0.0366 

+ 0.0010 

0.80 

1.20 

+ 0.0314 

- 1.2724 



+ 0.0760 

+ 0.0347 

0.80 

1.40 

+ 0.0801 

- 1.6374 



+ 0.1374 

. 

0.80 

- 0.02 

- 0.2981 

+ 1.1928 



- 0.2563 

- 0.3011 • 

1.00 

0.54 

- 0.2520 

- 0.0436 



- 0.2335 

- 0.2519 

1.00 

0.75 

- 0.1561 

- 0.4426 



- 0.1452 

- 0.1550 

1.00 

1.00 

- 0.0389 

- 0.9050 



- 0.0362 

- 0.0366 

1.00 

1.10 

- 0.0038 

- 1.0902 



- 0.0035 

- 0.0011 

1.00 

1.20 

+ 0.0015 

- 1.2752 



+ 0.0257 

+ 0.0047 

1.00 

1.35 

+ 0.0343 

- 1.5564 



+ 0.0621 

+ 0.0382 

1.00 

0.00 

- 0.3079 

+ 1.4728 



- 0.2306 

- 0.3116 

1.20 

0.58 

- 0.1515 

- 0.0826 



- 0.1502 

- 0.1513 

1 ..20 

0.80 

- 0.1022 

- 0.5116 



- 0.0945 

- 0.1009 

1.20 

1.05 

- 0.0426 

- 0.9890 



- 0.0277 

- 0,0401 

1.20 

1.15 

- 0.0214 

- 1.1792 



- 0.0037 

- 0.0184 

1.20 

1.22 

- 0.0080 

- 1.3126 



+ 0.0117 

- 0 . 0Q48 

1.20 

1.30 

+ 0.0059 

- 1.4650 



+ 0.0279 

+ 0.0096 

1.20 

1.40 

+ 0.0212 

- 1.6566 



+ 0.0462 

+ 0.0253 

1.20 

- 0.04 

/ - 0.2140 

+ 1.2460 



- 0 . 0401 ' 

- 0.0245 
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Table 
(k = -0.5 


Computation, of the profile in the 
generated by the function ( 8) . 
computed for points where M' = 0 


A 

7\ - oC 

9 

^v 

-0.194 

-0.294 

0.70 

0.3669 

-0.060 

-0.160 

0.40 

0.3986 

-1.000 

-1.100 

1.20 

0.1251 

-0.060 

-0.160 

0.30 

0.1044 

-0.094 

-0.194 

0.50 

0.2309 

-0.600 

-0.700 

1.00 

0.4385 

-1.410 

-1.510 

1.30 

-0.0165 

-2.000 

-2.100 

1.30 

-0.0902 

-1.000 

-1.100 

-0.03 

6.9885 

-1.700 

-1.800 

-0.10 

-0.0062 

-2.000 

-2.100 

-0.15 

-0.0792 

-0.720 

-0.820 

-0.02 

-12.9584 

-0.094 

-0.194 

0.36 

0.3748 

-0.720 

-0.820 

-0.02 

-12 . 9584 

-0.094 

-0.194 

0.36 

0.3748 

-0.197 

-0.297 

0.51 

0.5404 

-0.396 

-0.496 

0.78 

0.5105 

-1.000 

-1.100 

1.08 

0.1594 

-1.410 

-1.510 

1.17 

-0.0153 


0.000 


- 0.100 


0.00 
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III ' 

ij 4-0.1) 

physical plane. The flow is 


The values given below are 
(approximately) . 


.M's 

dx/dA 

dy/dA 

0.8476 

-0.5699 

-0.4800 

0.8505 

-0.4339 

-0.1834 

0 . 3829 

-0.3485 

-0.8962 

1.1894 

-0.5190 

-O'. 1605 

0.9525 

-O'. 4740 

-0.2590 

0.5050 

-0.6850 

-1.0669 

0.2453 

-0.2385 

-0.8590 

0.0272 

-0.5305 

-1.9110 

-1.8565 

64.3002 

-1.9298 

-0.1291 

0.6075 

-0.0609 

-0.0572 

1.0191 

-0.1540 

-146.8974 

243.4653 

-4.8703 

1.0478 

-0.5913 

-0*2259 

-146.8974 

243.4653 

-4.8703 

1.0478 

-0.5913 

| 

-0.2259 

0.8531 

-0.7775 

-0.4400 

0.6564 

-0.8152 1 

-0.8065 

0.2950 

-0.4009 

-0.7593 

0.1368 

-0.1885 

-0.4448 

0.6125 
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Figs. 1,2, 3, 4 




Figure 4.- Nomogram for the products (28). cos a = l/o.(M) . 

If_Oa = T(M,0i) , Ob = T(M,©2) , ...j then OA = 
Q(M)l(M,0i), OB = Q(M)T(M,©2) , ... .An obvious adjustment 
is mads when Q(M) < 1. By this procedure we may obtain any 
of the products q( v 'Tu. 
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Fig. 10 







Figure 12. 






